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Abstract 

We study direct CP violation of r ^ Ktt{i], r]')^ decays. By studying the forward and backward 
asymmetry, the interference of L = and L = 1 ampHtudes of the hadronic system can be extracted. 
By including the scalar and vector mesons into the chiral Lagrangian, we compute the form factors 
which correspond to L = and L = 1 angular momentum state of the hadronic system. We include 
real and imaginary parts of the one loop corrections to the self-energies of the scalar and vector 
mesons. The direct CP violation of the forward and backward asymmetry is computed using a 
two Higgs doublet model in which a new CP violating phase is introduced. We show how the CP 
violation of the forward and the backward asymmetry may depend on the new CP violating phase 
and the strong phase of the form factors. 

PACS numbers: 13.35.Dx.ll.30.Er,12.39.Fc,12.60.Fr 
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I. INTRODUCTION 



Two B factories, both Belle and Babar have accumulated the very large samples of r 
decays. The r lepton physics beyond the standard model, such as r lepton number violation 
and time reversal violation through the measurement of electric dipole moment have been 
studied. The CP violation of the hadronic r decay also has been investigated both theoreti- 
cally [1] and experimentally [2]. Various angular distributions including the quantities using 
the r spin polarization have been also proposed [l, ^. Recently, Belle and Babar reported 
the precise measurements of the branching fractions of r ^ KgTiv \}\ and K^tt^u sI. 

The improved measurement of the branching fraction for r — > Kriu has been also obtained 
jo]. About the r Ktxv decays, the hadronic invariant mass spectrum has been measured. 

Motivated by these measurements, we study the direct CP violation of K^Pi'{v) 
with P = 7r°, 1], 1]' and KsTt^uIu) . Non-vanishing direct CP violation in the processes 

may arise with some new source of the CP violation in addition to Kobayashi Maskawa 
phase 0] and the strong phase shifts of the final states of hadrons. 

In r KPu decays, the hadronic system KP may have the angular momentum L = 0(s 
wave) and L = l(p wave). The interference term of them can be extracted from the forward 
and backward asymmetry jsj. In the present paper, we define the asymmetry as the difference 
of the numbers of events for K meson scattered into the forward and backward directions 
with respect to the incoming r momentum in the hadronic CM frame. By comparing the 
forward and backward asymmetries for CP conjugate processes, the direct CP violation can 
be defined. The s wave and p wave of the hadronic amplitudes are related to the scalar 
and vector form factors in the time like region which have their own strong phases. To 
evaluate them, we use a chiral Lagrangian including scalar and vector meson resonances 
such as /t(800) and -ft'* (890). We compute the both real and imaginary parts of the one 
loop corrections to the self-energy of the resonances and obtain the strong phase shifts. We 
include the pseudoscalar meson loop correction and scalar and vector meson loop correction. 
The latter may give an important contribution to the form factors at higher invariant mass 
regions above 1 GeV up to ~ 1.7(GeV). 

As a new physics effect, we study a two Higgs doublet model with non- minimal Yukawa 
couplings to the charged leptons. The two Higgs doublets contribute to the charged lepton 
mass through the Yukawa couplings. In the non-minimal model, the interaction of the 
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charged Higgs boson to the r lepton family can be CP violating. The interaction generates 
the amplitude tr VtlH^ i^TL{uLSR)s-wa.ve- The interference with the charged current 
interaction due to W boson exchanged diagram may lead to the direct CP violation which 
can be measured in the forward and backward asymmetry. 

The paper is organized as follows. In section II, we show the hadronic chiral Lagrangian 
including scalar and vector resonances. In section III, we derive the form factors. In section 
IV, by fixing the finite renormalization constants, we numerically evaluate the form factors 
and the hadroic invariant mass spectrum. In section V, we introduce the two Higgs dou- 
blet model and present the direct CP violation. Section VI is devoted to conclusion and 
discussion. 



II. CHIRAL LAGRANGIAN INCLUDING SCALAR AND VECTOR MESONS 

In this section, we show the chiral Lagrangian with vector and scalar resonances. The 
following aspects are the main feature of the chiral Lagrangian. 

• U(1)a breaking effect is taken into account so that we can apply the Lagrangian to r 
decays into the final states including Krj and Krj'. 

• SU(3) breaking of the vector mesons are taken into account. 

About the inclusion of the scalar resonances, we followed the approach of Ref. [q]. About 
the vector meson sector, our Lagrangian is equivalent to the one in Ref. [lo| except SU(3) 
breaking effect for vector mesons. The chiral Lagrangian is given by, 

C = ^TiDUDU^ + BTtM{U + U^) - i^2pTr (^M^ - ^^M^^) r/o - ^r^o 
+ TiDf^SD^S - M^TiS^ 

- ^TrF^.^F^'^ + M'yTiiV, - + g^vTiS{V, - ^f, (1) 
2 9 9 

where S and V are the scalar nonets and vector nonets respectively. (See appendix A.) U 

is the chiral field and is given as U = exp(2i7r//) = vr is SU(3) octet pseudo Nambu 

Goldstone boson and r/o corresponds to U(1)a pseudoscalar of which mass is denoted by Mq. 

The U(1)a symmetry is broken by the mass term explicitly. The covariant derivatives for 
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the chiral field and the scalar field are given as, 

D^U = {d^ + iAL,)U, (2a) 

D^S = d^S + i[a^,S\, (2b) 



,0 ^ i!M±iM!, (2d) 

where denotes the external vector field corresponds to 811(3)^. M in Eq.([T]) is the chiral 
breaking term for the light quarks and is given by, 

M = diag.{mu,md,ms) 

= m, ■diag.(A,Arf,l). (3) 

A denotes — and A^ = IM^ in this work, we work in the isospin limit, m„ = m^. Below 

ITls Tils 

we explain how we determine the parameters in the Lagrangian of Eq.([T]). 

In the isospin limit, the vacuum expectation values of the scalar fields are given as, 

^01 - '^02 - -TT^, ^03 - -TT^, [^) 

which leads to the SU(3) breaking of the wave function renoramalization constants 
and the decay constants of the pseudo Nambu Goldstone bosons, 

ry 1 , Soi + Soj , . 

% = 1 + ^1 2/2 ' 

Fk = /v^, = /v^. (5b) 
The decay constant for rjs is written as. 
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where R is The srenerahzed Gell-Mann Oakes Renner relation becomes, 

m^F| = {rriu + TUs) {2B + g2{Soi + S03)) , (7a) 

mlF^ = 2m„ (25 + 2^2^01) , (7b) 

which can be used to express g2, gi and B with Eq.(jl]) and Eq.([S]) in terms of the 
physical quantities as, 



mj^ ]_ m| vl - A 

Y i+A 2AR2 

1 1 + A , 1 A , . . 

77 and r/' mesons and octet and singlet mixing angle ^os 

g2p in Eq.([T]) leads to the r/o and r/g mixing. The mass matrix for r^o and r^g sector is 
diagonalized as, 

.o,^>,.)("^:''*)(*)^-i(.-/)(^^-? °J| ;|, (9a, 




708 



-sin6'og 1 


II 




cos 6'og J 







(9b) 



where Mgg = ^s^mu ^j^^ ^g^j^g ^j-^g convention Mpg < 0. Beacuse the octet 
mass Mgg is given by 

Ml, = ^ (l^^i^l - ^M'^F' + M'^F') ' (10) 

M^g can be determined by Fk, F^, MxMn, and A. With M|g given by Eq. 1^, the 
parameters Mqq and Mgg are also determined by using the masses of rj and t]' as, 

M^o = Ml + Ml,-Ml,, (11a) 
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Ml^ = -sjMl,Ml - M^M'^,. (lib) 

Therefore, one may predict the ?7o and 773 mixing angle with the relation, 

1 2|M^ I 

The prediction of the mixing angle is rather close to the one experimentally extracted 
from J/ip ^ iviv') decays. (See for example {lll|.) 

^08th = -22.38(A = ^), -21.49(A = ^), (13a) 



^oscxp = -arctan^/ 1 —2 



^ -(22.36lii?)^ (13b) 
Vector meson mass spectrum 

The vector meson masses are given by the following formulae. 

M'y^,=M^y+g,y^^^. (14) 

With this formulae, the U(3) nonets vector mesons masses are given by, 

= Ml = Ml + g^vS,i, (15a) 

r2 



= M^ + (7iy5o3, (15b) 

A #2 ;\/r2 I "^01 + 'S'oa , s 

M^. = M^+^iy . (15c) 



One can fix the parameter giy as, 

91V = 2 (16) 
where A^* is the difference of the vaccumm expectation values in Eq.(jl]), 

Ag = 5„3-goi- "^'"tr^' . (17) 
One can also derive the following relation by using Eq. (fT5|) . 



M, = ^2Mi. - M| (18) 

The relation leads to the prediction Mp = 743 MeV which is about —4% smaller than 
the measured value. In table [U we summarize the numerical values for the parameters 
in the chiral Lagrangian of Eq. ([1]) . 
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= (MeV) 


800 


840 


760 


52 (MeV^) 


2.65 X 10^ 


2.79 X 10^ 


2.52 X lO'^ 


ffi(MeV) 


215 


225 


204 


Bnis (MeV^) 


9.24 X 10^ 


9.24 X 10^ 


9.24 X 10*^ 


TYls 


1 


1 


1 


AS = So3 - SoiiUeV) 


39.8 


37.9 


41.9 


5iv(MeV) 


12200 


12800 


11600 


9 


5.90 


5.90 


5.90 



TABLE I: The numerical values for the parameters in the chiral Lagrangian. We use Mk = 
494(MeV), = 135(MeV), Fx = 113 (MeV) = 92.2(MeV) and Tk- = 50.8(MeV) as input. 
g is determined with the width of K* . 



III. FORM FACTORS 

The hadronic form factors relevant for the processes r"' 



uK+P{P = 7r°,r7,V) are, 



{K+{pk)P{pp)\uj,s\0) 



F^'^mq'^ + ( Fr^'iQ') - ^F^'^iQ^) ] (19) 



with Qf^ = {px + Pp)^ and A^p = "n^x — rn'j,. The form factor denoted by F is the vector 
form factor and Fg is the scalar form factor. The form factors have been computed by using 

□ nnn 

the variety of the methods, Ref. [121. 1131 . 1141 . 1151 ]. In this work, we have used the hadronic 
chiral Lagrangian including the vector and the scalar resonances in Eq.([T]). We compute the 
loop corrections to the self-energy of the vector and the scalar resonances. The real part of 
the self-energy is divergent and we need to subtract the divergence. Corresponding to the 
subtractions, we have added the polynomials. Some of the coefficients of the polynomials 
are determined by the pole positions and the residues of the propagator for the resonances. 

To compute the form factors, let us write the V-A charged current in terms of hadrons. 
By differentiating Eq.([T]) with the external vector fields A^, we obtain the current as. 



t{^[s,d,s]e 



(20) 
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FIG. 1: The Feynman diagrams contributing to the form factors 



We first show the results of the form factors for Ktt final state. 

1 r R + R-^ (A5)2 



V2 



+ 



2FkF^ 2g^FKF^ 



+ 



2g^FKF^ 



'1 



A 



R 



1 - 



K* 



A 



R 



(21a) 



+ 
+ 

+ 
+ 



1 [r-^-R+ ^kABrDr-CI) 



2^2 I 



q^FkF^ Ar [Ar + Q^Br)Dr - Q^Cl 



1 



D, 



-9kK-k{Q'^) 



AS 



Ar+Q^Br 



{Ar + Q^Br)Dr - Q^Cl 
2M*2 



K ^ . 2Ax7r^L 

i : J-i 



VS 



(21b) 



2x/2F^F^^72 V Ar ' ■ Ar 

The form factors include the contribution of the Feynman diagrams shown in Fig. [H In 
Fig. [H the propagators for K* and k mesons are represented by the thick solid lines which 
include the one loop corrections to the self-energy. Let us consider the propagators for K* 
and K. They are obtained by inverting the inverse propagators for K* and k. 



{9'"'Ar{Q')+Q>^Q-^Br{Q^)) Q>^Cr{Q^) 
Q^Cr{Q^) Dr{Q^) 



)l 










I n ) 







(22) 
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where and J are source terms for K* and k respectively. Ajig'^^ + BrQ^Q'^ is the inverse 
propagator for K* and Dr denotes the inverse propagator for the /t. Cr denotes the mixing 
between the K* and k. Inverting Eq. (122]) . one can obtain the propagator. 



1 Kf"] 









g^-A^- + "q2 ^{An+Q-'BR)Dj^-Q-^Cj, {Ar+Q-^Br)Dr-Q'^CI | ^ 



^^lu A-l , Db^ 



(AH+Q2Sfl)DH-Q2C^ 



(Ar+Q2Br)Dr-Q2cI 



(23) 

To obtain the contributions to the form factors from Feynman diagrams Fig.(l-a) ~ Fig.(l- 
e), we set i = s, j = u in Eq. (120|) . 



ULltxSL 



Zi5 



2\/2 

2V ^ ^i^x 2g^FKF^ ^ ^ ^F^ 



2g^FKF^ 



1 



(24) 



V2 

For the diagram in Fig.(l-a), the direct couphng of the charged current to Kt^ can be easily 
extracted from Eq. (12^ . For the other diagrams, the amplitudes are obtained by multiplying 
the propagators in Eq. fl23|) with the production amplitudes of K* and k and the amplitudes 
corresponding to their decays. 

The matrix elements corresponding to Fig.(l-a) ~ Fig.(l-e) are given as, 



|0) 
|0) 
|0) 

|0) 
|0) 



1 f R-^-R f R-' + R Ml-. + g\^Sf 



(1-b) 

(1-c) 

(1-d) 

(1-c) 



2V2g^FKF^ 
ASg^K.iQ' 



2 ''^ \ 2 

1 



R^R 



CI 



Ar {Ar + Q^Br)Dr-Q''C^ 



2ni ] ' 



2V2FkF^ ^ 



2V2gFKF^ (Ar + g^i?,^)^^ - Q^Cl^^' 

Ml, g.KAQ'')CR 

2V2gFKF^ (Ar + Q^Br)Dr - Q^C], 



(25) 



9 




V K* 




(2-a) 



(2-b) 



(2-c) 



FIG. 2: The Feynman diagrams of scalar and vector mesons loop whcili contribute to the form 
factor for (i^P|'U7'^s|0). They can be written in terms of the self-energy correction function liys- 

To derive Eq. ( l25l) . we have used the production amplitudes for K* and n due to the vector 
current M7^,s, 



{K;\u^,s\0) = g,,V2^ 
(/€|m7^s|0) = Q^V2AS. 



(26) 



We also have used the strong interaction vertices which are given as, 

{K+7r^\CsPp\K+) = 



with q = Pk — Ptt- gnKwiQ'^) is the strong coupling for k — Kn defined by, 
9kkAQ^) = 9i— TT g2iSmu + ms) + Ak^{AS). 



(27) 



(28) 



In addition to the pseudoscalar loops, we have taken into account the vector and scalar 
meson loops denoted by Ily^. Each contribution is given by. 



{K+7l'>\uYs\0)\(^2-a) 

(ir+7rO|M7'^s|0)|(2_6) 
{K+n^\uYs\0)\^2-c) 



vs 



2^/2FKF^g^' 



A 



Q pQu \ 1 



where Ilys is identical to the self-energy function in Fig. (l3]-d). 



n 



VS — yy Q2 )^^VS + Q2 ^^VS- 



(29) 



(30) 



By denoting the self-energy corrections as 6Ai{, 6Br, Cr and corresponding to the 
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FIG. 3: (3-a): The pseudoscalar meson loop corrections to the self-energy for K* {6Aii,6BFi,). 
(3-b): The self-energy for k (SDr) . (3-c): The mixing term {Cr). (3-d): n^"^. 

Feynman diagrams in FigJSl one obtains 

Ar{s) = Ml, -s-iMK*rK*{s) + '^e-{SAR{s)) + Uls + {ao + ko) + {ai + ki)s, 



+ i 







487rs2 


167rs2 J 






487rs2 





08 



8/ 

2 



) ^9l*K. (^) sin^ ^08 + sRe.{5Bn{s)) + - I^y, 



Cr{s) = Re.(Cij(s)) + Co + ft^'l ^gt,K7r{s)gK*Kn 

~' 167rg2 4F ^2 ^'^^'?(^)^^'^^cos ^08 

^ 8 

. ^Kri'{s) Akv'Ftt / n .2/1 

Dr{s) = s- Ml + iM,T,{s) + Re.{6DR{s)) + do + dis + d2s\ (31) 

where s = and gK*K-K = A^gp^'p • '^'^^ momentum dependent widthes Fx*(s) and Fk(s) 
are given as, 



167r sM« \AFkF^ 

167r sM, ""^ 167r sM, J \4VSFkFs 
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where vkp = a/s^ — 2s(m\ + rn?p) + A^p (P = tt, r/, r/') and is related to the momentum of 
kaon pk in the hadronic rest frame as px = The real part of the self-energy corrections 
are divergent. We have subtracted the divergences at zero momentum as, 

5Ar{s) = 6A{s) - 6A{0) - sA'{0), 

6Bn{s) = 6B{s) - 6B{0), 

Cr{s) = Cis) - C{0), 

6Dr{s) = 5D{s) - SD{0) - sD'{0) - —D"{Q). (33) 

Then we have added the polynomials with respect to s which coefficients are finite renor- 
malization constants. We have added the polynomial oq + (i\s which corresponds to the 
twice subtractions for bAji in Eq. (1331) . For bBji and Cr, we have added a finite constant to 
each denoted by 60 a^nd Cq respectively. About the self-energy correction of the scalar meson 
6Dji{s), we need to subtract divergences up to s^. Therefore, we need to add the polynomial 
do + dis + (^25^ which is quadratic with respect to s. 

The self-energy corrections in Eq. ( l33l) . 6Ar ~ Cr, are given by, 

SAnis) = 

I (-^kAI + log ) + ^(-1 + log ^— ) - 2sRk. 

- E,4 + log^) + %±^log4 + £(-l + log^) 



- 2sRk^ - s ^—3 '- log + sm Uosl — ] [-2sRKr 

- SKr,'(o+iog r) + Q -1 + iog r) + — jx ^log — 

- s — + .0.2 y^^K + 5m^' - 22m j,m,,. 




(34) 



12 



SBr{s 
^9 



K*Ktt 



167r2 



cos^ 9, 



08 



sin^ e. 



08 



— log 

3 IJlKlTlTr 



Si 



Ktt 



El. 



6AL 



log 



m 



K 



mil 



4 + m\mi + mi 
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+ 



-^log — - — 
3 mxmjj 

-J log — - — 
3 mxmrfi 



^Kr, . m\ A m\ + m|^m ^ + 
^ 6Ai^ ^ ml 9 Ai 



+ 



m2, 



(35) 




m|f log -f s 



2A 



cos 6', 



1 1 1 "ll- 

m|^m^ log 



■'Kr\ 



A2 



sin ^08g^^2^ ^2 ^ 



( m'i-ml, log ^ \ 



V 



Kri' 



2Ax^/ 



(36) 



where Rpq and i?pQ are defined as, 

RPQ = dx [x^ - X{1 + ^) + ^) log 

R'pQ = ^RpQ - + RpQ, 



(37) 



with 



R^^^ = J^' dx log - x{l + ^) + 



M2 



16 



(38) 



We give the explicit forms for i?pQ (n = ~ 2) in appendix B. The inverse propagator for 
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the scalar meson is, 

1 

T 



sin^ 9n 



s' (2g,KA0)g'^KA0)J'^v'i^) + lg.KAOfJL,,{0)^ I , (39) 



where 



^PQl , g((Mp + MQ)^-g)g(5-(Mp-MQ)^) f^ApQ SpQ M2 



2 



+2 + I tan-i^ + tan-'i±^ ] } , (40) 

with SpQ = ml + ni^Q and i/|>q = - 2sEpq + A|,q = (s - (Mp + Mq)2)(s - (Mp - Mq)^). 
We also note, 



The self-energy corrections due to the vector and the scalar meson loop are also divergent. 



9lV ri ,A-d / " ^ ^ mf, 



. {2nYi {k^ - ml){{Q - kf - ml) 

= -% E Cvsig^^'K^ - Q'Q^'K^), (42) 
vs 

where Cys are factors determined by scalar and vector mesons which contribute to the loop 
and are given as, 

113 3 
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We subtracted the divergences of Ki and K2 as, 



K^vs = Ms) - KM - K'Ms, 

K2VS = K2{s)-K2i0). (44) 
Using Kivs and K2vsy we write the self-energies Ily^ and Uyg in Eq. flSUi) as, 

^vs = Cvs^^ivs, 



vs 

-\-\L _ Y\T 2 



5^Cy5i^2y5. (45) 



vs 

The exphcit forms for Kivs {i = 1? 2) are given as, 

- /-^ _ g + Ayg , g p(2) _ , g , ^VS _ ^VS 

s f_3 2 M^ + M^Ml + Ml \ MyMs f,^. s 

12M2 V 2 ^ 3 Al^s J s V + 12M2 

, My / Ey5 + 2M2 M2 s + - 3M2 M^Sys 

+ log — — -r -TTT- H r — ^ + 



Ms \Avs 4Ay5M2 sAvs 12M2 



/I // ^ 



1 1 M^log^ 11 M|logi^ 
- 16^Af:;\ 3 18^^^ 3 + ~9 

My^t , ,,4 ,,2 ,.2,,2A..„1..^n 1 Z,(2) 



+ M^M| - M^M|Av5 log ^} - (46) 



2 

The absorptive parts of Kivs (^ = 1)2) are written as 



^vs 5 Ayg Syg ^ 

327rs^ 6M2 2M2 GM^ GsM^ ^' 

^yg /H _ ^yg 3Ay5 A^ 
487rsM2/ 2s 2s s^ 



I-(A-.v.) = - + ^ + (47) 



In Aji and of Eq. fl32|) . we have added the polynomial lis to 11^^ — Ily^. and ko + kis to 
Hys- "^^^ polynomials corresponds to the subtraction of the divergent parts of Eq. (jSj). 

The above derivation of the form factors for Ktt final state is easily extended to Kt] and 
Krj' case. 



v/2 I 2 2FkFs 

nT / 2M2 \ M2 / M2 _ , ^ ^ 



2f72F^F8 V An I 2g^FKFs \ A 
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Q 



1 (n^-,^^ 

?>gFKF^ {Ar + Q^Br)D - Q^C^ 



ttT / 0/1^2 



(49a) 



1 (n'^^^l_ 

Ar+q^Br 

g.K^'{Q^)Ml. - 3M|-.A^,.A^ Cr 

SqFkFs (Ar + Q'Br)Dr - Q^C\ 



with i?8 = and A^^(/) = m|. — ""^^(/j- To derive Eq. (148 ap and Eq. (I49al) . we used the 
strong interaction vertices for K* Kr] and K* — ^ Krj' , 

{K^r,\Cvpp\K:-') = cos^os^^fV^^- 

{K^V'\^VPP\K+) = smOos—^VSq., (50) 
where q = Pk — Pp {P = Vi v')- We also use k — » Krj^''^ vertices which are given by, 

a^KrjC) (Q^) = 91 ^ ^2 (5m, - m„) - SA^^^o (AS). (51) 
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rriK 


493.7 


'/ 


957.8 




135.0 




775.5 




891.7 


Mr 


1777 




1019 


Tr 


2.26 X 10-9 




782.7 




0.2257 




547.5 


Tk* 


50.8 



TABLE II: The numerical values for masses, widths and Vus used in our numerical analysis. The 
units of mass and width are MeV. 

IV. NUMERICAL ANALYSIS OF THE FORM FACTORS 



In this section, we give the numerical results of the form factors. We summarize the 
numerical values for the masses of hadrons and width which are used for the numerical 
analysis in table [Tll We first determine the finite renormalization constants. There are ten 
constants. We first renormalize Re.A/j and Re.Di? so that each inverse propagator has zero 
at the on-shell mass s = Mj^, and s = respectively. And then we require the residues 
of the propagators Re.^ and Re^^j^ are unity on their pole masses. Because of the four 
conditions, we can constrain the parameters ao + /cq, cti + ^i, c^O) di and We show the four 
conditions below. 



ao + A;o + (ai + ki)Ml. = -Re.Mi?(M|. 

T 

( ;, 1 1 1-: . I f ( A-1 o —I— I I 

tti + ki 



Re.nTg(M|0, 



dRe.{6An + Illsj 



do + diMl + d2Ml = -Re.SDn 
2 , (iRe.5DR 



di + 2M^d2 



ds 



=M2- 



(52) 



Moreover we set four parameters Co,bQ,li,d2 to be zeros. Below we show the numerical 
values of constants for = 800 MeV case. 



ao + ko 
ai + ki 
di 
do 



-5.68 X lO^(MeV^), 

0.345, 

-0.140, 

-1.87 X lO^(MeV^). 



(53) 
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(MeV) 


760 




800 


840 


oi + ki 


0.33961 




0.34471 


0.35100 


«0 + ^0 


-63519.7 




-56822.9 


-52534.2 


di 


-0.139154 




-0.140467 


-0.134301 


do 


6485.53 




-1871.23 


-16362.1 




-4.16 X 10^ 




-4.03 X 10^ 


-3.92 X 10^ 


h 


0.6558 




0.6561 


0.6564 


Br(K7r°) 


0.00416079 




0.00416079 


0.00415823 


Br(i^?7) 


0.000162063 




0.000162015 


0.000162325 


Br(Kr/') 


4.14285 X 10"^ 




3.87288 X 10^6 


3.94976 X 10-^ 



TABLE III: The numerical values for the finite renormalization constants. They are chosen so that 
the the branching fractions of r — > Kni/ and r K-qu can be reproduced. 

We also show the values of the constants in Table ilTTl for = 760 MeV and 840 MeV. Be- 
cause the other parameters were set to be zero, there are only two undetermined parameters. 
We choose ko and ki as the parameters to be adjusted. They are fixed so that the branching 
fractions for r i^vrz/ and r if?7z/ can be reporoduced. For {kQ, ki) = (—4.03x10^,0.656), 
we obtain 

5r(r^ ^ K^nu) = 0.416 x 10-^ 

Br{T^ K^r]u) = 1.62 x 10'^ (54) 

which are close to the expeimental results and [sl, 

Br{T- K'TT^u) = 0.416 ± 0.003 ± 0.018 x 10"^ 

Br{T- K-r]u) = (1.62 ± 0.05 ± 0.09) x 10~^ (55) 

The branching fraction for r — > Kri'u becomes, 

Br(r Kri'u) = 3.87 x 10"^ (56) 

We have plotted the hadronic invariant mass spectrum for Rtt^Kt] and Kt]' cases in Fig. HI 
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The formulae can be found in 

((^ + l)PK'\F'''is)\^ + ^\Fs^'is)\^) , (57) 

where px is the momentum of kaon in the hadronic CM frame. In the hadronic invariant 
mass spectrum for K^n^ at low invariant mass region, K* resonance can be seen. Just below 
K*, we can see the effect of k(800). At the high invariant mass region, the new thresholds 
due to the vector and scalar channels are open and these effects can be seen in Kit, Kt] and 
Krj' cases in Fig. HI 

dBr 



. 00001 




1 . X 10"^ 




1 5 » 


1 . X 10"^ 




I ^ 1 
1 ■ ^ 

! » 
■ ! * 
1 1 1 




1 . X 10"^^ 







600 800 1000 1200 1400 1600 1800 

Vs(MeV) 



FIG. 4: The hadronic invariant mass spectrum for Ktt^ (solid line), Kr] (dashed line) and Krj' 
(short dashed line) cases. We choose = 800 MeV and the other parameters are fixed as in the 
corresponding columns of Table Hill 



V. FORWARD AND BACKWARD ASYMMETRY AND CP VIOLATION 



In this section, we compute the forward and backward asymmetry [8[ and the direct CP 
violation for r — KPu decay. We start with the double differential rate of the unpolarized 
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r decay 



dBr 



1 G^IK/K 



dy/sd COS 6* 



m 



^cos^^ + sin^^ 



^KP\2 



m 



(58) 

where Q is the scattering angle of kaon with respect to the incoming r in the hadronic CM 
frame. The forward and backward asymmetry extracts the interference term of the vector 
form factor and the scalar form factor. 

[,'rfcosg ^/f , - fVcosg- 

J dJsd cos ti J ~\ I 



dy/sd c 



dBr 



Pk 



^ \FKP\ COS 6^^ 



3 / m2 2 I ^ I 



y 3s 



(59) 



with 6^^ = arg.(wF). As we can see from Eq. (IS^ . the forward and the backward asymme- 

^ s 

try is determined by the ratio of the scalar and the vector form factors. It is also proportional 
to cosine of the strong phase shift 6st- In Fig. [5l we show the vector and the scalar form 
factors and their ratio for Kit case. We also show the strong phase shifts in Fig. [61 The 
forward and backward asymmetries for Kit , Krj and Krj' cases are shown in Fig. [71 As can 
be seen from Fig. [71 the forward and backward asymmetry for Kit case is large near the 
threshold region where the scalar contribution is also large. (See Fig. [5l) We can expect 
about 50% asymmetry for r —>■ Knu decay. 




600 800 1000 1200 1400 1600 1800 



00 300 1000 1200 1400 1600 1800 



Vs(MeV) 



Vs(MeV) 



FIG. 5: Left: The vector form factor (thick solid line) and the scalar form factor 10 x \F, 

(thin solid line). Right: The ratio j^fer 



KttI 



r 



By including new physics source of CP violation , we can predict the direct CP violation 
of the forward and backward asymmetry. To be definite, we start with non-minimal two 
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600 800 1000 1200 1400 1600 1800 

Vs(MeV) 



FIG. 6: The phase of the vector form factor: = aig.F (short dashed hne) and the phase 
of the scalar form factor 6^'" = avg.F^'^ (dashed hne). The strong phase shift 6^^'" = 6y'^ — 6^'" 
is shown with thick sohd hne. and are shown with long short dashed line and solid line 
respectively. 
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FIG. 7: The predictions of the forward and backward asymmtries of r — > Kiriy (thick solid line) 
r — > K'qiy (dashed line) and r Krj'u (short dashed line) in the standard model. 

Higgs doublet model in which the two Higgs doublets have couplings to the charge leptons, 
— ^ = yiijemHi Lij + y2ije^HlLLj + y^jdmHi Qij 

+ y2i3URiH2QLj + y2ijl^iH2LLj + ^^V^i{vRi)'' + h.C, (60) 

where we assume that H2 is coupled with neutrinos and both Hi and H2 are coupled with 
the charged leptons. The neutrino mass is given by the seesaw mechanism. However, the 
right-handed Majorana neutrino with the mass M^v much heavier than the electro- weak 
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scale does not affect on the interaction terms with the mass dimension equal to four at all. 
Therefore, we keep the terms which are not suppressed by a factor in discussion below. 
In two Higgs doublet model, without loss of generality, one can parametrize the Higgs fields 



H2 = e' ^ \ \ Hi = 1T2H: = e * 2 , (61) 

where 6c p is the relative phase of the vaccuum expectaion values of Higgs and its value can 
be determined from Higgs potential. The charged current interactions of the model are. 



cos (3 V2Mw 



gmdi tan p _ 
V2Mw 



- -^v:,W-duruL - -^V-Zl^JuW^^ + Ml^W+^^W., + h.c, (62) 

where the charged lepton masses mi = diagonal(me, m^, m,-) are obtained by the diagonal- 
ization. 

VR^{-y,v,e''-^ + y2V2e-''-^)Vl = mi. (63) 

Eq. ( |63ll can be used to express yi in terms of the charged lepton mass and the other Yukawa 
coupling j/2. By introducing, 

Y2 = Vny2Vl (64) 

one can obtain Eq. ( 162|) . The four fermi interactions induced by the charged Higgs exchange 
are, 

= ^^{^(^--^^-ii^)"^'^-^' + ^^^4 

X {dV^M {^d{l - 75) + cot^ (3{1 + 75)) m| + h.c, (65) 

where r2ij denotes the non-minimal couplings of charged Higgs boson between the charged 
lepton Ij to the neutrino z/,, 

^2ij = ^ySM > (66) 
i 
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where Y^^ denotes the standard model Yukawa couphngs for charged leptons, ^f^^ 
niij. The amphtude of the two Higgs doublet model is, 



Amp.(r- ^ U,K-P) = -^v:^u,r{l - l,)u^ 

[SU^,-^Q,)F +(Ml-^tan/5) + ^— ^r..jF, Q 
where the matrix element of the scalar current is given by, 



(67) 



{K-Asu\0) = -^Fr. (68) 
iiig 1 1 III 

Then r KPvi branching fraction is, 

— d^sd cose — - f 2%^^;^ 1 \— ^ + V ^^'^ ' + t'^-- ' 



2 



- 5,.^p^cos^^Re.(F^^F,^r)|, (69) 
where we define 

- {^,r (l - gta.^ ,) + (70) 

for i = e, fi and r. We neglect the small corrections proportional to up quark mass. For CP 
conjugate processes K^n^Ui is obrained by replacing r2 in the amplitude r~ K^tt^Ui 

with its complex conjugate Therefore the direct CP violation of the forward and backward 
asymmetries is given as, 

A A '^fs^^i^^^^st ^Q^tan^/jy 

where we parametrize CP violating phase of the flavor diagonal coupling as, 

r2rr=\r2rr\e'''-^. (72) 

We set the flavor off-diagonal couplings in r2ir to be zeros. We note that in the isospin 
limit m„ = m^, the CP asymmetry of Eq. 071 p for K^tt^ case is identical to the direct CP 
violation of r~ — Kg'K'Vi and r"*" Ksii^Vi. Contrary to the CP violation of the total 



branching ratios which is sensitive to the CP violation of the mixing of and 16| . 
the CP violation of the forward and backward asymmetries does not depend on the mixing 
induced CP violation. 
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The CP violation of the forward and backward asymmetry is shown for K^ir^ case in 
Fig. [HI By taking \r2TT\ = 1, one can see that CP asymmetries can be as large as a few 
%. At low invariant mass region y/s < 900 (Me V), the direct CP violation is negative while 
at high invariant mass region ^/s > 900 (MeV) , the CP violation is positive. The sign 
is correlated to sin^^j as can be seen in Eq. ( 17T1) . From Fig. [6], we can see sin^^t also 
changes its sign around y/s = 900 (MeV). We also change the charged Higgs boson mass. 
For Mh > 500 (GeV), CP violation is suppressed to less than 1 %. In Fig. [9], we also show 
the CP asymmetries for Krj and Ki]' cases. We note the sign of the CP violation is opposite 
to the sign at high invariant mass region of Kn case. 



0.08 
0.06 
0.04 
0.02 
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-0 . 04 
-0.06 
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FIG. 8: CP violation for the forward and backward asymmetries of r ^ Kttv. The charged Higgs 
boson mass is changed as MniGeY) = 200 (thick solid line), 250 (dashed line), 500 (short dashed 
Une) and 700 (sohd Une). The other parameters are tan/? = 50, |r2Tr| = 1 and 92tt = f • 

VI. CONCLUSION AND DISCUSSION 

We have studied CP violation of r — > KPv {P = Tx^^rj^rj') decays and r — >■ KsTtu. CP 
violation for the forward and backward asymmetries is computed using the two Higgs doublet 
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FIG. 9: CP violation for the forward and backward asymmetries of r ^ Kni' (thick sohd hne), 
r — > K'qv (dashed Une) and r Kr]'v (short dashed hne). We choose the parameters as Mh = 300, 
tan/3 = 50 \r2TT\ = 1, ^2rT = f 

model which has the non-minimal Yukawa couplings to the charged leptons. The effect of 
the CP violation is numerically studied. The direct CP violation depends on the strong 
phase shifts sin(5v' — ^s)- To evaluate the phase shifts, we use the chiral lagrangian including 
the vector and the scalar resonances and take into account of the one loop corrections. We 
compute both real part and imaginary part of their self-energies. The divergences of the real 
part is subtracted properly. We have determined the finite renormalization constants so that 
the branching fractions for Kti and Kr] modes are reproduced. We also take into account 
of U(1)a breaking so that r^o and rj^ mixing can be incorporated. With those improved 
treatment, we can predict the hadronic invariant mass spectrum and CP violation even at 
high invariant mass region for Kt[ case as well as the same observables of Krj and Krj' cases. 

The effect on the direct CP violation of the non-minimal coupling is studied. For the 
non-minimal Yukawa coupling as large as that of the standard model Yukawa coupling of 
r lepton, we have predicted a few % CP asymmetries within the parameter region with 
Mh = 200 ~ 300(GeV) and with tan/3 ~ 50. 
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APPENDIX A: SCALAR AND VECTOR MESON NONETS 



Here we show the scalar and vector meson nonets. 
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APPENDIX B: FUNCTIONS H 



The function i?pg is defined as, 

Jo s s 



in) 
PQ 



(Al) 



(Bl) 



There are two regions for s of interests depending on below the threshold, i.e., (1) {Mp 

Mqf <s< {Mp + Mqf or the above threshold (2) (Mp + Mqf > s. 
For the case (2), 

-1 



(n) 
PQ 



a;"log((a; — a){x — h) — ie)dx, 



(B2) 
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with a and b are given as, 



2s ' 



2s ' 

with < a, 6 < 1. We show the real part of R^pg, 
Re-inPg) = -2 + log(l-a)(l-6)-alog^^-Mog^, 



a 

' ^ - — log h - logfl - a) \ log + - 

2 ^l-a 2^^ ' 4 2 2 ^1-6 2 



/ Ti(2) \ 6a^ + 3a + 2 1 , x , a 

Re.(i?>A) = — + -log(l-a) + — log- 



^PQ, - -3-6^^ 



662 + 36 + 2 1 63 ^ 

+ -log(l -6) + — log- 



For the case (1), 



18 3 °' ' 3 °l-6 



R^pI = [ dxx''log{{x-(3f + a^) 
Jo 



' ^ dy{P + yriog{y' + a'), 



where, 



a 



^{{Mp + MQY-s){{s-{Mp-MQf) ^ \I'^Iq 
2s 2s 

s + ApQ 



/3 = 



2s 

We define the indefinite integrals. 
Using the integrals, one can write, 

R?Q - rPgil -(5)- rPQ{-(5) + (5R%, 



The indefinite integrals are given by, 

^PQ = \iy^ log(y2 + o;2) - ^ + 2a^y - 2o? arctan ^), 
^3 3 a 

= \ + log(y^ + c?) - y') , 

rpQ — y log(y2 + a^) — 2y + 2q; arctan 
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